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Abstract 

The type Hi unprojection is, by definition, the generic complete in- 
tersection type II unprojection, in the sense of [P] Section 3.1, for the 
parameter value k = 1, and depends on a parameter n > 2. Our main 
results are the explicit calculation of the linear relations of the type Hi 
unprojection for any value n > 2 (Theorem 3.16) and the explicit cal- 
culation of the quadratic equation for the case n — 3 (Theorem 4.1). In 
addition, Section 5 contains applications to algebraic geometry, while 
Section 6 contains the Macaulay 2 code for the type Hi unprojection 
for the parameter value n = 3. 

I Introduction 

The first appearance of the type II unprojection was in the study of elliptic 
involutions between Fano 3- fold hypersurfaces in [CPR], while the theoret- 
ical foundations were developed in [P] using valuations. Valuations were 
useful to prove the existence of certain relations, but didn't provide any ex- 
plicit formulas for them. The present paper is an effort towards the explicit 
calculations for the type II unprojection. 

We define the type II& unprojection to be the generic complete intersec- 
tion type II unprojection, in the sense of [P] Section 3.1, for the parameter 
value k > 1. According to [P], it depends on a parameter n > 2, increases 
the codimension from nk — 1 to nk — 1 + {k + 1) and preserves Gorensteiness. 
To our knowledge, the only previously done explicit calculation for the type 

II unprojection was for the type Hi case for n = 2 ([CPR], [R]), which is 
reproduced for completeness in Subsection 4.4 below. 

Section 3 contains the calculation, using homological and multilinear 
algebra, of the linear equations of the type Hi unprojection for any value 
n > 2. The main result is Theorem 3.16, which provides explicit formulas 
for the linear relations. 
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Section 4 contains the calculation of the quadratic equation of the type 
Hi unprojection for the case n = 3. The main result is Theorem 4.1 which 
provides an explicit symmetric - in the sense of SL3 invariance, see Sub- 
section 4.2 - formula for the quadratic relation for this case. The proof of 
Theorem 4.1 is based on the explicit equality (4.6) which was verified using 
the computer algebra program Macaulay 2 [GS93-08]. In Subsection 4.3 we 
briefly sketch how we arrived to the formulas contained in Theorem 4.1 by 
computer assisted calculations (using the computer algebra program Maple) . 

As an application of the above results, we sketch in Section 5 the con- 
struction of two codimension 4 Fano 3-folds. The quasismoothness checking, 
using the explicit equations obtained in Sections 3 and 4, was done by the 
computer algebra program Singular [GPS01]. 

The calculation of the quadratic equation for the type Hi unprojection 
for n > 4 remains open. We believe that the explicit formulas of the lin- 
ear relations obtained in the present work together with representation and 
invariant theoretic techniques could, perhaps, lead to their calculation, and 
also to a better understanding of the complicated looking formula (4.2). 
To our knowledge, the problem of calculating the type 11^ unprojection for 
k > 2 is also open. 

ACKNOWLEDGEMENTS: I wish to thank Janko Boehm, Stephen Donkin, 
Sebastien Jansou, Nondas Kechagias, Frank-Olaf Schreyer and Bart Van 
Steirteghem for useful discussions. Parts of this work were financially sup- 
ported by the Deutschen Forschungsgemeinschaft Schr 307/4-2 and by the 
Portuguese Fundagao para a Ciencia e a Tecnologia through Grant 
SFRH/BPD/22846/2005. 

2 Notation 

As already mentioned in Section 1 by type Hi unprojection we mean the 
generic complete intersection type II unprojection in the sense of [P] Sec- 
tion 3.1 with fixed parameter value k = 1. According to [P], the type Hi 
unprojection depends on an integral parameter n > 2, the initial data for 
the unprojection is specified by a triple 

Ix C Id C O amb (2.1) 

(where Ix and Id are ideals of O am b), and the unprojection constructs an 
ideal 

Iy cO amb [s , Sl ], (2.2) 
where sq,si are new variables. 
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Fix n > 2. The ambient ring O am b will be the polynomial ring 

Oamb = ^[xi, ■ ■ ■ , x n , j/i, . . . , y n , z, A P j,Bf m ] 

where 1 < p < n — 1, 1 < i < j < n and 1 < I < m < n. The ring O a m6 
corresponds to the ring denoted by 9 amb in [P] Section 3.1, the variable Xj 
corresponds to the variable aij in [P] Section 2, and the variable yi corre- 
sponds to the variable a>2i in [P] Section 2. (In Section 4 O am b will change 
slightly to allow 2 to be invertible.) 
We denote by M the 2 x 2n matrix 

M= A/i ... //,. z Xl ... zx n \ 
\xi ... x n yi ... y n J 

which corresponds to the matrix with the same name in [P] Equation (2.1). 

The ideal Id C O am b in (2.1) is the ideal generated by the 2x2 minors 
of M and corresponds to the ideal I g D in [P] Section 3.1. 

For I < p < n — 1, we denote by A p the n x n skew-symmetric matrix 
with (ij)-entry (for i < j) equal to A\- (and zero diagonal entries), and by 
B p the n x n symmetric matrix with (ij)-entry (for i < j) equal to Bfj. 

For 1 < p < n — 1 we set 

n 

l<i<j<n i=l l<i<j<n 

where, 

-Fjj = Xiyj — Xjyi, Gij = lix/iyj zx{Xj), G mm = y m zx m 
for 1 < i < j < n and 1 < m < n. In matrix notation 

f p = xAPy* + yBPy 1 - zxBPxK (2.4) 
The ideal Ix in (2.1) is the ideal 

Ix = (/l, • • • , fn-l) 

and corresponds to the ideal I 9 X in [P] Section 3.1. 
In addition, we set 

x = (xi, . . . ,x n ), y = (yi, . . . ,y n ). (2.5) 
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Remark 2.1 According to [P] Proposition 2.16, the ideal ly of O am b[so, si] 
contains Ix as a subset and is generated by Ix together with n polynomials 
[sq, si] affine linear in s , S\ of the form 

k = zxiS + ViSi + a}, 

where 1 < i < n and aj G O am b, together with n polynomials l n +i, ■ ■ ■ , hn £ 
Cambf^O) s i] affine linear in so, si of the form 

k+n = ViSo + XiSi + of, 

where 1 < i < n and <r| G O am i,, together with a single affine quadratic 
polynomial q G Oam6[so; «i] of the form 

(/ = — zSq + lower terms, 

where 'lower terms' means an affine linear polynomial in S0) s i; i- e -> °f t ne 
form eis + e 2 si + e 3 with ei, e 2 , e 3 G O amb . 
In other words, we have the equality 

I Y = (fl,...,fn-l) + (h,...,l2n) + (q) (2.6) 

of ideals of O am b[so, s\]. Calculating ly means providing explicit formulas 
for h,h, ■ ■ ■ ,hn and q. 

For simplicity, we call the elements of Ix the original relations, we call 
h,. .. ,hn the linear relations and, finally, we call q the quadratic relation. 

For the parameter value n = 2 the ideal ly has been explicitly calculated 
by Reid in [R] Section 9.5. For completeness, we reproduce his results in 
Subsection 4.4. 

In Theorem 3.16 we calculate, for any n > 2, the linear relations h,. . ■ , hn- 
In addition, in Theorem 4.1 we calculate the quadratic polynomial q for the 
parameter value n = 3. 

The calculation of the quadratic polynomial q for n > 4 remains open. 

3 Linear relations of Type Hi unprojection 

In this section we use the notations of Section 2. 
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3.1 Generalities 

Assume R is a ring, L is an i?-module, and h: L — > R a homomorphism of 
i?-modules. For p > 1, [BH] p. 43 defines an i?-homomorphism d p h : A P L — > 
A^L by 

p 

cf /t(ui A • • • A Up) = ^(-f ) i+1 /i(uj)ui A • • • A Uj A • • • A Up, (3.1) 
i=i 

for all ui, . . . , v p G L. (For p = 1 we set A°L = i? and d^h = h.) The maps 
d p h define the Koszul complex 

. . . ^ A P L a p-i l _ . . . ^ A 2 L ^ L R o 
associated to the homomorphism h. 

Proposition 3.1 Assume hi,h 2 - L — > R are two homomorphisms of R- 
modules. For p > 2 we have the equality 

a?- 1 ^ o d p h 2 + a?- 1 h 2 o <Fhi = 0, 

of maps A P L A p - 2 L. 

Proof Indeed, 

= dP-\h! + h 2 ) o d p {hi + h 2 ) = df- x h x o d p h 2 + o cfpi, 

since 

(F^/ll O (F/ll = Cf _1 /l 2 O d P h 2 = 0. 

QED 

3.2 The second complex 

Assume L is a free O am ;,-module of rank n, and let ei, . . . , e n be a fixed basis 
of L. We define four am b-homomorphisms h±, . . . , : L — > am b, with 

hi(ei) = yi, h 2 (ei) = x { , h 3 (ei) = zx u h 4 (ei) = y u 

for 1 < i < n. In addition, we define, for 1 < p < n, homomorphisms 

<t> p : A p L® A P L -> A P ~ 1 L A P ~ 1 L 
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(where by definition A°L = O am b) as follows: 

Ma,b) = Q = ((Fh^a) +dPh 3 (b),dPh 2 (a) +dPh 4 (b)), 

for a,b e A P L. 

Proposition 3.2 For p > 2 we have 

4> P -i ° 4>p = o. 

Proof Using the equalities /14 = hi and /13 = zh 2 the result follows from 
Proposition 3.1. QED 

The proof of the following proposition will be given in the Subsection 3.3. 

Proposition 3.3 The complex 

L* : O amb © O amb L © L A 2 L © A 2 L A 3 L © A 3 L <- . . . (3.2) 
is exact. 

3.3 Proof of Proposition 3.3 

The proof of Proposition 3.3 will be based on the Buchsbaum-Eisenbud 
acyclicity criterion as stated in [BH] Theorem 1.4.13. We first need the 
following combinatorial lemma. 

Lemma 3.4 Let p be an integer with 1 < p < n. Then 




Proof We have 

n , n. n-1 , <. 

^(-iy-v( n ) = (-ir* + x;(-ir p ( •) 

^-ir+EVinglJj + f: 1 )) 
= (-ir- + (;: 1 1 ) + (-ir l - = (;: 1 1 ), 
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where we used the Pascal's rule (™) = + ( n j X ) which is valid whenever 
1 < j < n - 1. QED 

For the following we fix an integer p with 1 < p < n. We set 

n 

r p = ^(-l)^rank (A P L A P L), 
j=p 

where for a finitely generated free O am b-module N we denote by rank N the 
minimal number of generators of N as O am b-module. Using Lemma 3.4 we 
get 

0.3) 

Denote by I rp (<fip) the r p -th Fitting ideal of the map (fi p , that is the ideal 
of O am b generated by the r p x r p minors of any matrix representation of 
4> p . (For more information about Fitting ideals see, for example, [BH] Sec- 
tion 1.4.) 

We will need the following general property of Koszul complexes. 

Lemma 3.5 Fix t G {1,2} and s,p with 1 < s,p < n. Denote by Mf the 
matrix representation of the map 

d p h t : A p L ^ A P ~ 1 L 

with respect to the bases A • • • A ei p of A P L and A • • • A ei p _ 1 of A P ~ 1 L 
induced by the basis (ei,...,e n ) of L. For a suitable ordering of the two 
bases there exists an (r p /2) x (r p /2) submatrix of Mf which is diagonal with 
diagonal entries equal to at or —at, where a\ = y s and 02 = x s . 

Proof Consider the subset of the above mentioned basis of A P L consisting 
of the elements where e s appears on the wedge. This subset has 
elements, hence (r p /2) elements using (3.3). The result follows by restricting 
the homomorphism d p ht to the O am fe-submodule of A P L generated by this 
subset and noticing that the corresponding matrix has a diagonal submatrix 
of the desired form. QED 

Proposition 3.6 Fix s,p with with 1 < s,p < n. There exists a polynomial 
g s G I rp {4>p) of the form 

g s = y7+gl (3.4) 

where g\ involves only the variables y s , z,x s and has degree in y s strictly less 
than r v . 
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Proof It follows immediately from the definition of 4> p and Lemma 3.5. 
QED 

Proposition 3.7 Let F be either the field Q of rational numbers or a finite 
field Ti/(q), where q is a prime. The ideal I of O am b ®7L F generated by the 
image of I rp {4>p) under the natural map O am b — > O am b(g>iF has codimension 
greater or equal than n. 

Proof Using Proposition 3.6 it is clear that there exists a monomial order- 
ing for the variables yj, z,Xj (with indices 1 < j < n) such that the leading 
term of the polynomial g s appearing in (3.4) is equal to y s p , for 1 < s < n. 
Since over the field F an ideal / and the leading term ideal of / have the 
same codimension the result follows. QED 

Recall (cf. [BH] Section 1.2) that the grade of a proper ideal / C O am b 
is the length of a maximal O am fe-regular sequence contained in /. 

Proposition 3.8 The ideal L rp (4> p ) has grade greater or equal than n. 

Proof Combining Proposition 3.7 and [L] Theorem 3.12 we get that I Tp (0 P ) 
has codimension in O am b greater or equal than n. Being a polynomial ring 
over the integers, the ring O am b is Cohen-Macaulay. Using [BH] Corol- 
lary 2.14 the result follows. QED 

We now finish the proof of Proposition 3.3. Combining Lemma 3.4 
and Proposition 3.8, Proposition 3.3 follows from the Buchsbaum-Eisenbud 
acyclicity criterion as stated in [BH] Theorem 1.4.13. 

3.4 The first complex 

Fix n > 2. Let TV be a free O am fe-module of rank n — 1, with basis 

e l j • • • ) e n— 1 • 

We define an C am b-homomorphism ip : N — > O am b, with tp(e p ) = f p , for 
1 < p < n — 1, where f p was defined in (2.4). As in Subsection 3.1, we have 
the Koszul complex 

N * : 0amb ^NS^A 2 N<^A 3 N^... (3.5) 

Proposition 3.9 The sequence f±, . . . , fn—i is an (D am \,-regular sequence. 
As a consequence, the complex N* is exact. 
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Proof For an element g of O amo we denote by r)(g) the result of substitut- 
ing to g zero for z, zero for A^- with 1 < p < n — 1 and 1 < i < j 1 < n, and 
zero for B^- for 1 < p < n — 1 and 1 < i < j < n with ^ (p,p)- That 
is, we set zero the variable z and all possible A\- and all possible Bfj with 
the exception of B pp (for 1 < p < n — 1). 
It is clear that for 1 < p < n — 1 we have 

V(fp) = B p pp (y p ) 2 . 

Therefore, the ideal (rj(fi),...,r)(f n -i)) C O am b has codimension in O am b 
equal to n — 1. As a consequence, the ideal (/i, . . . , f n -i) of O am f, has also 
codimension n—1. Since O a mfe is Cohen-Macaulay the sequence /i, . . . , / n -i 
is an O am 5-regular sequence. The exactness of the Koszul complex N* fol- 
lows from [BH] Corollary 1.6.14. QED 



3.5 The first commutative square 

We define two O am fe-homomorphisms ui,U2- N — > L, that will make the 
following diagram 



Oamb 



Oamb © @amb *~ 



iV 



- L © L 

commutative. That is, they will satisfy the two conditions 



1p = hi O U\ + z/l2 ° ^2 
= hi o U2 + tl2 o Ui 



We set, for 1 < j < m, 



(3.6) 
(3.7) 



and u 2 (ej) = ^c^e,, 



with 



and 



n i( e i) = E c ^ e * 

i=l i=l 
j— 1 rain 

^ = - E + E 4** + E B ^ + E 4^ 

Z=l l=i+l 1=1 l=i+l 

i n 

4 = - E B l*i - E 

1=1 l=i+l 
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for 1 < i < n and 1 < j < n — 1. 

In matrix notation, u\ and U2 are given, with respect to the above defined 
bases (ej) of N and (e^) of L, by the following two nx(n-l) matrices: 

Ul = (-AV + B 1 y t | - AV + B 2 y t \ ...\- A n ~ 1 x t + B^y 1 ) 

u 2 = (-B 1 x t | - B 2 x t | ... | - B'^x 1 ). (3.8) 

That is, for 1 < j < n — 1, ui is the matrix with j'-th column equal to 
— A J x* + B^y 1 , while 112 is the matrix with j-th column equal to —B^x 1 . 

Proposition 3.10 The maps u\ and U2 satisfy the equations (3.6) and 
(3.7). 

Proof An easy calculation using (3.8). QED 
3.6 The basic formula 

The key result of the present subsection is Proposition 3.12, essentially a 
formal consequence of (3.7), which relates the differentials and the alter- 
nating products of the homomorphisms hi and Uj. We need the following 
definitions. 

Assume ai, ct2 : N — >■ L are two am ft-homomorphisms. 

For p>0we have a natural induced homomorphism A p a\ : A P N — ► A P L. 
By definition, A°ai is the identity map O am b — > O am b, For p > 1, the map 
A p a\ is uniquely specified by the property 

A p ai(c\ A • • • A Cp) = ai(ci) A • • • A a\(c p ) 

for all ci, . . . ,Cp G AT. We denote the map A p ai also by alt(a^). 
We will now define for p, q > 1 an am b-homomorphism 

alt«,a«): A p + 9 N -> A P+9 L. 

We set 

alt (a^a^) = u>/, 

where the sum is over all subsets / C {1, . . . , p + (/}, with |/| = p and, by 
definition, wi = wi A ■ ■ ■ A w p+q with Wi = a± if i G I and w i = 02 if i ^ 
More precisely, even though the tensor product map 
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does not induce a map A p+q N — > A p+q L, taking the sum of those maps 
over the family of subsets / as above induces a well-defined map A p+q N — > 
A p+q L, and this is the map denoted by alt(af , a|). 

Remark 3.11 For p = q = 1 we have 

alt(ai,a 2 )(ci A c 2 ) = ai(ci) A a 2 (c 2 ) + a 2 (ci) A a 2 (c 2 ), 

for all ci, c 2 G iV. We will also denote alt(ai, a 2 ) by ai A a 2 + a 2 A a±. 
For p = 2, q = 1 we have 

alt(cii, a 2 ) = ai A ai A a 2 + ai A a 2 A ai + o 2 A ai A ai, 

in the sense that 

alt(af, a 2 )(ci A c 2 A c 3 ) = ai(ci) A ai(c 2 ) A a 2 (c 3 ) + 

ai(ci) A a 2 (c 2 ) A ai(c 3 ) + a 2 (ci) A ai(c 2 ) A ai(c 3 ) 

for all ci, c 2 , c 3 £ iV. 

For the following proposition, which will be proved in Subsection 3.7, 
recall that ui, t/ 2 : N — > L were defined in Subsection 3.5, and hi, hi : L — > 
Camb were defined in Subsection 3.2. We also use the notational conventions 
alt (u^, itSj ) = alt(u^), alt (1/5,1*2) = alt^), alt(tij )_2 , 1X3) = for p = 1 and 
q > 1, and alt(n^, «2~ ) = for q = 1 and p > 1. 

Proposition 3.12 a) Let u: M — > L and /i: L — > O am b fee too (arbitrary) 
homomorphisms of O am \,-modules. Assume p > 2. PFe Ziawe 

d p - l h o alt^" 1 ) = alt(n p - 2 ) o d p -\h o u). 

b ) For p,q > 1 we /iai>e 

(d p+ ^ 1 /i 2 ) oaltK,^" 1 ) + (dP+'-^i) oaltK" 1 ,^) = (3.9) 
alt«~ 2 , u q 2 ) o d p+q - 1 (h 1 o m) + alt«, uf 2 ) o o u 2 ), 

where the equality is, of course, as maps A p+q ~ 1 N — > A p+q ~ 2 L. 

3.7 Proof of Proposition 3.12 

We will use the following two lemmas. 
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Lemma 3.13 Assume p, q > 1 and u±,U2- N — > L are two O am b-module 
homomorphisms. We have 

alt(n^, u q ) = ui A alt(n^ _1 , u|) + u 2 A alt^, nrj! *)) 

in i/ie sense f/iaf 

alt(uf , n^Xc A ci) = 

ui(c) A alt(n^, M2)(ci) + n 2 (c)A alt(n^, n2 _1 )(ci) 

for all c £ N and a £ A**' -1 ^ 

Proof Immediate from the definitions. QED 



Lemma 3.14 Assume p,q > 1, ui,u 2 - N — > L are two O am b-module ho- 
momorphisms, and h: N — > O am b is an O arn t,-homomorphism. We have 

alt« u\) o (d p+<?+1 /i) = h A alt« ul) - m A [alt«~\ u«) o 

-u 2 A [alt«,^ _1 ) ort 

in the sense that 
alt«, u«) o (d p+q+1 h)] (c A ci) = 

h(c) [alt«, u|)(ci)l - ui(c) A (alt« _1 , u«) o d p+q h){a) 
-u 2 {c) A (alt«, o a? +q h){ Cl ) 

for all c £ N and a £ A p+q N. 
Proof Since 

0? +q+1 h(c A a) = h(c) Cl - c A d p+q h( Cl ), (3.10) 

we have, using Lemma 3.13, that 

alt«, u\) o (d p+q+1 h)(c A ci) = /i(c)(alt«, n^)(ci)) 

- ui A alt« _1 , u«) + u 2 A alt«, u« _1 )l (c A d p+q h( Cl )) 

and the result follows. QED 
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We will now prove part a) of Proposition 3.12. Assume p > 2. For 
ci, . . . , Cp-i G N we have 

o alt(u p_1 )(ci A • • • A Cp_i) = d p_1 /i(u(ci) A • • • A u(c p _i)) 
p-i 

= ^2(- 1 ) i ~ lh «d))(u(c 1 ) A • • • A u(ci) • • • A n(c p _i)) 
i=i 

= alt(u p ~ 2 ) o dP" 1 ^ o u)(ci A • • • A Cp_i). 

We will now prove part b) of Proposition 3.12. Assume that p, q > 1, 
that AT, L are O am b-modules and that ui,u 2 : N — > L and h\,h 2 : L — > O am b 
are four O am (,-module homomorphisms with the property 

h\ o U2 + h,2 ° u\ = 0, (3-11) 

cf. (3.7). We will show by induction onp + g that (3.9) holds. If p = g = 1 
the result is clear, since it is exactly our hypothesis (3.11). 

Assume now that p + q > 3 and that (3.9) holds for the (p — 1, q) and 
(p, q — 1) cases. 

Using Lemma 3.13 we have 

dP +q - 1 h 2 o altK,^ -1 ) + dP+^hi o alt« -1 ,^) = 

cP +q - 1 h 2 o [«! A alt^" 1 , ^r 1 ) + u 2 A alt«, t^" 2 )] 
o [ Ul A alt«~ 2 , u q 2 ) + u 2 A alt« _1 , uJf 1 )] 

which, using (3.10), is equal to 

(h 2 o ui) A ait(u?~\ ul' 1 ) — ui A [d p+q - 2 h 2 o altK -1 ,^" 1 )] 
+ (h 2 o u 2 ) A alt«, u q 2 ' 2 ) - u 2 A [(F+^ 2 /i 2 o alt«, u^f 2 )] 
+ (/ii o ui) A alt«~ 2 , tig) - ui A [<F +9_2 /n o alt(nr 2 , A)\ 
+{hi o « 2 ) A ait(u?~\ ul' 1 ) - u 2 A [d p+9 - 2 /i 1 o altK -1 ,^" 1 )] 

which, using (3.11), is equal to 

-ui A [d p+q - 2 h 2 o altK -1 ,^" 1 ) + d p+q ' 2 hi o alt«~ 2 ,^)] 
-«2 A [d p+q - 2 h 2 o alt«, u^ 2 ) + d p+q - 2 hi o alt« _1 , ul" 1 )] 
+(/i2 o u 2 ) A alt(n^, «2~ 2 ) + (^i ° ^i) A alt(n^ 2 , u q ) 
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which, using the inductive hypothesis, is equal to 

(h 2 o u 2 ) A alt«,^- 2 ) - ui A [alt«-\^- 2 ) o d p+q - 2 (h 2 o u 2 )\ 

-U2 A [alt«, u* -3 ) o d p+q ~ 2 {h 2 o U2 )] 
+ o m) A alt«" 2 , u q ) - m A [alt« -3 , u«) o cf +9_2 (/ii o Ul )] 

-«2 A [alt(ur 2 ,«2 _1 ) ° d p+9_2 (/ii o m)] 
which, using Lemma 3.14, is equal to 

alt«" 2 ,«l) o o m) + alt«,u^ 2 ) o d p+q -\h 2 o U2 ) 

which finishes the proof of Proposition 3.12. 

3.8 The connecting homomorphisms 

For an integer p, with 1 < p < n — 1, we set 

/ alt(n^) + z alt(u^~ 2 , u 2 .) + z 2 alt(ti^ -4 , tx 4 ) + . . • \ 
p= \alt{ur\u 2 ) + z8lt(u{-\ul) + z 2 a\t{u p - 5 ,4) + .. .) ' 

with the summation as long as all exponents are nonnegative. More pre- 
cisely, using the notational convention z° = 1, and denoting, for j = 1, 2, by 
Zp the j-th row of T p , we have 

[f] [^] 
Tp = a \t(u p - 2 \ uf), T 2 =Y j z i a lt(u p - 2i -\ uf +1 ). 



i=0 i=0 



We consider T p as an O am ;,-homomorphism A P N — > A P L © A P L. Taking 
the two rows of T p we get two am (,-homomorphisms T p ,T 2 : A p N — > A P L. 
For example, 

uA rp _ ( A 2 u\+zA 2 u 2 \ _ f A 3 m + zalt(ui,u?,) 



'A 4 m + 2 alt(ti 2 , u|) + z 2 A 4 u 2 
■ alt(ti 3 , U2) + z alt (tii, it 3 ) 



T 5 = 



A 5 ui + z alt(n 3 , tt 2 ,) + z 2 alt (tti, u\) 
alt(tt 4 , ti 2 ) + z alt(ttf , ti 3 ) + z 2 A 5 u 2 J ' 



For the following proposition recall that the map <p p was defined in Sub- 
section 3.2 while the map ip was defined in Subsection 3.4. 
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Proposition 3.15 For any p, with 2 < p < n — 1, we have a commutative 
diagram 

dP4> 



A p - l N 
Tp-i 



A p_1 L © A^L 



A P N 
A P L © A P L 



That is, 

<f> p oT p = T p _i o d?il> 
as maps A P N -> A P ~ 1 L © A P_1 L. 

Proof Fix p > 2. For i G Z we set 

Ci = alt(ui,u?P) 

if < z < p and Q = otherwise, and we set 

Di = alt(«i,«5" 1_i ) 

if < z < p — 1 and Di = otherwise. 

For the rest of the proof all sums are for ieZ. 
We have 

By Proposition 3.12, for any igZwe have 

cf /ii o Ci_i + d p h 2 od = A-2 o <F(/ii o m) + A o d p (h 2 o u 2 ). 
Using (3.6), we need to show that 

d p h 1 o + z(F/i 2 o T p 2 = TpLi o cFQn o ui + z/i 2 o u 2 ) 

and that 

(F/i 2 o Tl + o T p 2 = T p 2 _! o <F(/ii o ui + zh 2 o u 2 ). 
Using (3.12) we get 

d p h 1 tf + zd p h 2 T* = ^[/(d p /iiC p _ 2l + d?h 2 C p - 2i+1 )] 



(3-12) 
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= Y^[z i (D p - 2i - 1 dP(h 1 u 1 ) + D p _- 2i+1 dP(h 2 u 2 ))} 

= Y,[^(D (p _ 1) _ 2i (d p (h 1 u 1 ) + zcP{h 2 u 2 )))} = Tl_ l( P{h lUl + zh 2 u 2 ). 
Similarly, 

<Ph 2 Tl + dPh^ = Yy^bvCp-n + efhiCp-M-i)] 
= ^ i {D p . 2i - 2 <F(h 1 u 1 ) + D p _ 2i d?(h 2 u 2 ))} 

= Y.lz'iDw-K-ii'Ffcui) + zd?(h 2 u 2 )))} = rf^cFQuu! + zh 2 u 2 ). 
QED 

3.9 Linear equations of type Hi unprojection 

Assume n > 2 is given. We use the notations of the previous subsections. In 
Subsection 3.8 we defined two O am fe-module homomorphisms T^_ x : 

A n ~ 1 N — > A n_1 L. 

There exist, for j = 1, 2 and 1 < z < n, unique element u\ E O a mb such 
that 

n 

I^-i(ei A • • • A e n _i) = ^(-l)*" 1 ^ d A ■ • ■ A ej A ■ ■ ■ A e n . 

i=i 

For 1 < i < n we define the polynomials li,k+ n £ Camfe^ s i] with 

Zj+n = yjSO + XiSi + of (3.13) 

Using Proposition 3.15 and arguing as in [R] Section 9.5 we get the 
following Theorem. 

Theorem 3.16 Fix a parameter value n > 2. The polynomials h, ■ ■ ■ ,l 2n G 
Oamb[so, si] specified in (3.13) are the linear relations of the type II\ unpro- 
jection in the sense of Remark 2.1. 
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4 Quadratic relation of Type Hi for n = 3 



Assume now that we are in the type Hi unprojection with n = 3. We 
follow the notations of Section 2 with only one change: Since our symmetric 
formula (4.2) for the quadratic relation will need 1/2 as coefficient, our 
ambient ring O am 5 will now be 



Oamb = Z[-][a;i,.. .,x 3 ,yi,.. . ,y 3 ,z,A^,Bf m ], 



with indices 1 < p < 2 and 1 < i < j > < 3, 1 < / < m < 3. 

Recall that in Section 2 we defined two symmetric 3x3 matrices B 1 and 
B 2 by 

(b(i BJ i2 b{ 3 \ 

BJ = B{ 2 B{ 3 , 

\sym B 3 33 J 

for j = 1, 2, two skew-symmetric 3x3 matrices A 1 and A 2 by 

(0 A\ 2 A\ 3 
A{ 3 
-sym 

for j = 1,2, and two 1x3 matrices x, y with 

x = {x 1 ,x 2 ,x 3 ), y = {yi,y 2 ,y 3 )- 

We now define three induced 3x3 symmetric matrices ad B l , &dB 2 and 
adS 1 * 2 . 

For a square nxn matrix M with entries in a commutative ring i? with 
unit, we denote by ad M the nxn matrix with fe/-th entry, for 1 < k, I < n, 
equal to (—l) k+l times the determinant of the submatrix of M obtained by 
deleting the Z-th. row and the k-th. column. It is well-known that we have 
the identities M(adM) = (ad M)M = (detM)/ n , where I n is the identity 
nxn matrix over R. 

In our situation, since B 1 and B 2 are symmetric 3x3 matrices we have 
that ad-B 1 and ad-B 2 are symmetric 3x3 matrices, such that for j = 1,2 

f Bi 2 B j 33 -{Bi 3 ) 2 -{B{ 2 Bi 3 -B{ 3 B{ 3 ) b{ 2 b{ 3 - b{ 3 B{ 2 \ 
bAB>=\ Bi ± Bi 3 - (B{ 3 ) 2 ~{B{iB 3 23 - B{ 3 B{ 2 ) . 

sym B{ 1 B 3 22 -{B{ 2 ) 2 ) 
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Notice that for 1 < k < 3, 1 < / < 3, we have 



nl nl \ 
J -'in J - > mj ) 



where {i,k,m} = {j, l,n} = {1,2,3} and i < m, j < n, and similarly for 
adB 2 . 

We define adB 1 ' 2 to be the 3x3 symmetric matrix, such that, for 1 < 
k < 3, 1 < I < 3, 

(adB^ 2 ) kl = (-l) k+l (BhBl n + BlBl„ - BLbL - BLBl 



where {i, k, m} = {j, I, n} 
definition we have 



{1,2,3} and i < m, j < n. By the above 



ad B 
adB 
adB 
ad B. 



1,2 

n 
1,2 
12 
1,2 
13 
1,2 
22 
1,2 



pi d2 i d2 pi 
i3 22 £ *33 ' i3 22 i3 33 



2^23^23 



-B\ 2 B 33 



d2 pi i pi p2 i p2 pi 
-D-19-D33 + ti \3 t >23 + ^13^23 



B 12^23 + B \2 B 23 



12-^33 
2 pi 



B 13 B 22 



'13 J 
B 13 B 22 



ad B 23 
adB 33 



nl d2 1 d2 pi o pi p2 

B nB 33 + B n B 33 - 2B 13 B 13 

pi p2 / >2 pi 1 pi p2 I p2 pi 

-B\\ B 23 ~ B n B 23 + B 13 B 12 + B 13 B 12 

pi p2 1 p2 pi o pi p2 
-011-022 + -011-022 ~~ 2 -Ol2-Oi2 

and the rest of the entries of adi? 1 ' 2 are determined by symmetry. 
We now define two polynomials a, b G O a m6- We set 





( Xi 


X2 


X3 


a = det 


^23 


-A 1 
^13 


A\2 




V^23 


-Ai 3 


A±2 



x i{A\ 2 A\ 3 — ^2^13) + X 2(A\ 2 A 23 — ^412^23) 



+ %3 (^13^23 ~~ ^13^2 



23^ 



and 



b = - 



yB 1 (adB 2 )B l y t + yB 2 (ad B 1 )B 2 y t 

+ 2xA l (ad B 2 )(A 1 ) t x t + 2xA 2 (ad B^A 2 )^ 1 

- 2x A 1 (ad B 1 ' 2 )(A 2 ) t x t 

- zxB 1 (ad B 2 )B l x t - zxB 2 {ad B l )B 2 x t 

+ AyB l (adB 2 )(A l ) t x t + AyB 2 (adB 1 )(A 2 ) t x t 

3 

- ( E 4(^% -^))( E 4 ad 4) 

l<i<j<3 l<i,j<3 

- ( A U x iyj- x jyM E B fj adB h) 

l<i<j<3 l<i,i<3 



(4.1) 



(4.2) 
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The two polynomials / 



i> h corresponding to the ones in (2.4), are 



fj = xA j y l + yB'y* - zxB j x l 



(4.3) 



for j = 1,2. 

We denote by u\ the 3x2 matrix 



Ul = (—A 1 x t + B 1 ^ 



A 2 x t + B 2 y t ) 



and by U2 the 3x2 matrix 



u 2 = {-B 1 x t 




These matrices correspond to the maps with the same name defined in (3.8). 
For j = 1,2, we set <^i(u J ) to be the 3x1 matrix 



The matrix (f>i{u^) corresponds to the map /\ 2 Uj and the matrix 4>2{u l ,u 2 ) 
corresponds to the map alt(ui,U2) which were defined in Subsection 3.6. 
We denote, for 1 < i < 3, 



where we use the usual notation Mij for the ij-th entry of a matrix M. 
Clearly, for 1 < i < 6, k £ O am b[so, si] is affine linear with respect to so 
and s\. The polynomials correspond to the polynomials with the 

same name in (3.13), hence by Theorem 3.16 they are the linear relations of 
the unprojection in the sense of Remark 2.1. 

The main result of this section is the following theorem, which specifies 
the quadratic relation of the unprojection ring for the parameter value n = 3. 
It will be proved in Subsection 4.1. 




and we also define the 3x1 matrix 




k = zxiS + ViSi + [4>i{u v ) + z<j)i{u 2 )]n 
-3 = ViSo + XiSi + [foiu 1 ,u 2 )]n, 



(4.4) 
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Theorem 4.1 The polynomial 

q = s j- zsq- as + b £ O amb [s , s ± ], (4.5) 

where a was defined in (4-1) and b was defined in (4-2), is the quadratic 
relation of the type 11\ unprojection for the parameter value n = 3 in the 
sense of Remark 2.1. 

In other words, the ideal ly C O am b[so, s±] defining, in the sense of 
Section 2, the type \\\ unprojection is equal to 

lY = (fi,h) + (h,...,h) + (q), 

where the two polynomials f\, fi generating Ix were defined in (4-3) and the 
six polynomials l±,. . . ,Iq were calculated in (4-4)- 

4.1 Proof of Theorem 4.1 

We first prove the following reduction lemma. 

Lemma 4.2 To prove Theorem 4-1 it is enough to prove that that the ele- 
ment (xi + %2 + x?)q of O am b[so, si] is inside the ideal (/i, /2, h, ■ . . , Iq) of 

Oamb [SO, Si]- 

Proof Indeed, by [P] Proposition 2.16 there exists q € O am b[so, si] of the 
form 

q = s l — zs + u 

where u is affine linear in so, si, such that the ideal of the unprojection ring 
is equal to J, with 

J = (h, f2,h, ■ ■ ■ , k, q) c O amb [s , si]. 

Since by the assumptions of the lemma (xi + X2 + x%)q G J, and by [P] Sec- 
tion 3.1 J is a prime ideal of O am h[so, si] we get q G J. Therefore, the afhne 
linear with respect to so, si element q — q is in J, as a consequence (compare 
proof of [P] Proposition 2.13) we get 

(fi,f2,h,---,k), 

hence 

J = (fi,f2,h,---,k,q) 
which finishes the proof of Lemma 4.2. QED 
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For the rest of the proof we define some more notation: 

We define two polynomials §{,§2, such that, for j £ {1,2}, we have 

gj = ( Xl B 3 n + x 2 B{ 2 + x 3 B{ 3 )(ad B p u + ad B\ 2 + ad B\ 3 ) 
+{xiB 3 21 + x 2 B\ 2 + x 3 ^ 3 )(adSf 1 + ad-B^ + ad5f 3 ) 
+(x 1 B J 31 + x 2 B 3 32 + x 3j B^ 3 )(ad5 3 J 1 + ad5f 2 + ad5f 3 ) 

where p£ {1,2} is uniquely specified by {j,p} = {1, 2}. 
Moreover, we define two 1x3 matrices 

L\ = (h, l 2 ,h), L2 = {h,hih) 

and two 3x3 matrices Mi and M 2 by 

( M i)ij = ^3(^1, B} 2 , B} 3 , Bj x , B] 2 , Bj 3 , B n ,B? 2 , Bf 3 , B^,B^ 2 , B 2 j3 ) 

and 

(M 2 )ij = 03 (A}i , A} 2 , A} 3 , A\y , A\ 2 , ^4^3 , Bj! , S| 2 , #J 3 , Bj! , B 2 j2 , B 2 j3 ) 
- <f>3 (Aa , A\ 2 , A? 3 , ^ , A 2 i2 , ^3 , B) x ,Bj 2 , B l j3 , ,Bj 2 , B) 3 ) 

where 

<M a i' a 2,a 3 ,&i,& 2 ,6 3 ,ci,c 2 ,c 3 ,di,d 2 ,d 3 ) = (ai(eZ 3 - d 2 ) + a 2 (di - d 3 ) 
+a 3 {d 2 - di)) - (ci(6 3 - 6 2 ) + c 2 (6i - 63) + c 3 (6 2 - &i))- 

Finally, we define the polynomial 

w = xM x (L x ) 1 + yMi(L 2 )* + zM 2 (L 2 )'. 

Proposition 4.3 We Ziaue i/te following equality 

2( Xl +x 2 + x 3 )q = -w + g 1 f 1 +g 2 f 2 -2s Q (l 1 + l 2 + l 3 ) + 2s 1 (l 4 + l 5 + k) (4.6) 

of elements of O amb [s , si]. 

Proof The verification was done using the computer algebra program 
Macaulay 2 [GS93-08]. QED 

The proof of Theorem 4.1 follows now by combining Proposition 4.3 with 
Lemma 4.2. 



21 



4.2 Invariance under SL 3 

Let P G SLa(Z) be a 3 x 3 matrix with integer coefficients and determinant 
one. Define two 1x3 matrices x = (xi,X2,Xs) and y = (2/1,2/2,3/3) by 

X = xP~ l , y = yP~ t } 

(where P~ l is, of course, the transpose matrix of the inverse of P), two 
3x3 skew-symmetric matrices A> and two 3x3 symmetric matrices P>i , 
for j = 1,2, by _ _ 

A? = P* A J P, BJ = P'B J P 

and finally set z = z. Moreover, we denote by Aj the Zm-th entry of the 

matrix A> and by Bj m the Zm-th entry of the matrix Bi . 

It can be checked that if we denote by a and b the polynomials of Expres- 
sions (4.1) and (4.2) respectively, with the variables without tilde replaced 
by the corresponding expressions with tilde, we have a = a and b = b. 

4.3 How we got to Theorem 4.1 

In this subsection we briefly sketch the method that allowed us to calculate, 
using the computer program Maple, the expressions for a and b appearing 
in Theorem 4.1. 

Step 1. Recall the explicit expressions for k in (4.4). We started from 
the expression 

2 [s (h + h + h) + si(U + h + k)] 

which is equal to 

2(xi +x 2 + x 3 )(sj - zsl) + 2H, 

where 

3 3 
H = -soJ^iCu 1 ) + zMu 2 )]n + si^2[Mu\u 2 )]n. 
i=l i=l 

It is easy to see that 2H is a homogeneous quadratic polynomial with re- 
spect to the variables A\- , . More precisely, 2H has a (unique) natural 
representation as a sum of terms of the form ^ij^ki^kl P ms sum 
of terms of the form -^\j^ti c t^i P ms sum °f terms of the form 

A lj B ll c ?jkl P lus sum of terms of the form B ij B kl c Fjkl- 

Step 2. We first handle the 'subpart' of 2H which is sum of terms of the 

form AjjA^cfj^. These term appear only as contributions from So(/i+/2+^3) 

and, moreover, this subpart is equal to 2aso, where a was defined in (4.1). 
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Step 3. We now handle the 'subpart' of 2H which is sum of terms of the 
form AjjB^Cy^. We notice that there is a kind of complementarity of terms 
for the pairs (Zi , Z4) , (h, h), (h, h) m the sense of the following example: 

A partial sum of the coefficient of A\ 3 Bf 2 coming from l\ and I4 is 
2yiX\SQ + 2x\s\ which can be written as 2xi{y\SQ + Using I4 we 

get an expression not involving sq and s±. Similarly, another partial sum 
coming from I2 and ^5 is 2x\ (2/2^0 + X2S1) and we can use £5 to to get an 
expression not involving so and s±. A similar procedure can be done for the 
the partial sum coming from l 3 and l%. 

By symmetry, the handling of the 'subpart' of 2H which is sum of terms 
of the form AfjB^cf^ is similar. 

The handling of the 'subpart' of 2H which is sum of terms of the form 
B}jBli c fjki ls similar, but we need to take care of preserving the symmetries. 
This is the only part where we actually need the coefficient 2 in 2H. 

To give an example, a partial sum of the coefficient of B^B^ coming 
from l\ and I4 is 

2 [-s (-xix 3 z - ym) + si(y 3 xi + yix 3 )} 

which we first write as 

2/3(2/150 + x 1 s 1 ) + xi(zx 3 s + y 3 si) + yi(j/3.so + + x 3 (zx 1 s + y 1 si) 

and then we use the appropriate linear equations U to get an expression not 
involving s and s±. 

Step 4- At the end of step 3 we arrived to an expression not involving 
so,si- By suitable subtraction (in a symmetric way) of multiples of fi and 
f2 we get a polynomial divisible by x\ + X2 + x 3 . The quotient is b, which 
after some further term by term effort can be written in the form (4.2). 

4.4 Type II4 unprojection for n = 2 

This subsection contains the explicit equations obtained by Miles Reid for 
the type Hi unprojection with parameter value n = 2. It is taken from [R] 
Section 9.5. 

For n = 2 the ambient ring is 

Oamb = %[xi, x 2 , 2/1, 2/2, z, A u , B n , B 12 , B 22 }, 
the matrix M corresponding to the one in (2.3) is equal to 

M = hi 2/2 zxi zx 2 \ 
\xi x 2 2/1 2/2 / ' 
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there is a unique polynomial 

/ = A 12 (x!y 2 - x 2 yi) + B u (yl - zx\) + 2B 12 (yiy 2 - zx\x 2 ) + £22(2/2 - zx 2 2 ) 

corresponding to (2.4), the linear equations h, . . . , I4 are given by 

h = zxi-sq + yisi + {xiA 12 + yiB 12 + y 2 B 22 ) 

l 2 = zx 2 s + y 2 si + (x 2 A 12 - yiBu - y 2 B 12 ) 

h = Viso + xisi + (-xiBi 2 - x 2 B 22 ) 

k = V2So + x 2 si + (xiBu + x 2 B 12 ) 

and the quadratic equation q is given by 

q = sj - zsl - A 12 s - {B l2 f + B U B 22 . 

In other words, using the notations of Section 2 we have the equality of 
ideals of O amb [s , si] 

lY = (f) + (h,...,k) + (q). 

Remark 4.4 Since / = y 2 l\ — y\l 2 + zx 2 k — zx\k, we even get 

I Y = (h,---,k) + (q)- 

Remark 4.5 By [P] the ideal Iy is Gorenstein codimension 3. It is easy to 
see that it is equal to the ideal generated by the 5 submaximal Pfaffians of 
the 5x5 skew-symmetric matrix (with entries in O am b[so, s\]) 

/0 xi x 2 yi y 2 \ 

-s -B 22 si + B 12 

-S1 + B12 -B11 

-zs - A 12 

\ -sym / 

(For a discussion about the Pfaffians of a skew-symmetric matrix see, for 
example, [BH] Section 3.4.) 

5 Applications to algebraic geometry 

We believe that the explicit formulas of Theorem 4.1 can be used together 
with Gavin Brown's online database of graded rings [Br] for the proof of the 
existence of a number of (singular) Fano 3-folds with anticanonical ring of 
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codimension 4. We discuss below two such examples, the first of which is 
due to Selma Altmok, and has the interesting property | — Kx\ = 0. We also 
expect that the explicit formulas of Theorem 4.1 together with the orbifold 
Riemann-Roch theorem obtained in [BS] can lead to the construction of 
new codimension 4 Calabi-Yau 3-folds. 

5.1 The example of Altinok 

This example is taken from [R] Example 9.14 and is due to Altinok. Consider 
the weighted projective space P 2 (l,3, 5) with coordinates u,v,w and the 
weighted projective space P 5 (2, 3, 4, 5, 6, 7) with coordinates x,v,y,w, z,t. 
We define the map P(l, 3, 5) -» P(2, 3, 4, 5, 6, 7) given by 

2 7 

x = u , v = v, y = uv, w = w, z = uw, t = u . 

The equations of the image D of the map are 

y z t xv xw x 4 \ ^ 
v w x 3 y z t J ~ ' 

that is, 

yw = vz, yx 3 = vt, zx 3 = wt, y 2 = v 2 x, yz = vwx 
yt = vx A , z 2 = w 2 x, zt = wx A , t 2 = x 7 . 

A general complete intersection X\2,i4, containing D is obtained by choos- 
ing two general combinations, one of degree 12 and the other of degree 14, 
of the above equations of D. 

We perform a type Hi unprojection of the pair D C X 12; i4 to get a 
codimension 4 3-fold 

Y C P 7 (2,3,4,5,6,7,8,9). 

Notice that the new variable so has degree 8 and the new variable s\ has 
degree 9. 

After substituting to the formulas (3.13) and (4.5) we checked the qua- 
sismoothness of Y using the computer algebra program Singular [GPS01]. 

5.2 The second Fano example 

In this subsection we sketch a construction, suggested by Brown's online 
database of graded rings [Br], of a codimension 4 Fano 3-fold 

Y CP 7 (1,1,2,2,2,2,3,3) 
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starting from a (nongeneric) co dimension 2 complete intersection Fano 3-fold 
14^^(1,1,2,2,2,3). 

Write x\, X2,yi,y2,U3, w for the coordinates of P 5 (l, 1, 2, 2, 2, 3), and con- 
sider the subscheme D C P 5 (l, 1, 2, 2, 2, 3) with equations 

rank ( Vl 2/2 W mXl 2/3X2 ^] < 1 
\xi x 2 y 3 yi y 2 w J 

By definition, the ideal Id of D is generated by the 2x2 minors of the 
above matrix. Denote by X±q C P 5 (l, 1, 2, 2, 2, 3) a general codimension 2 
complete intersection (with equations of degrees 4 and 6) containing D. The 
equations of X^q are obtained by choosing two general combinations, one 
of degree 4 and the other of degree 6, of the above equations of D. 

We perform a type Hi unprojection of the pair D C to get a codi- 
mension 4 3-fold Y C P 5 (l, 1, 2, 2, 2, 2, 3, 3). Notice that here the new vari- 
able so has degree 2 and the new variable s\ has degree 3. 

After substituting to the explicit formulas (3.13) and (4.5) we checked 
the quasismoothness of Y using the computer algebra program Singular 
[GPS01]. 

6 Appendix: Macaulay 2 code 

This section contains the Macaulay 2 code for the type Hi unprojection with 
n = 3. 

— M2 code for type II_1 unprojection for n=3 
kk = QQ 

S = kk [Alpl2 ) Alpl3,Alp23,A2pl2 ) A2pl3 ) A2p23 ) Blpll , Blpl2 , Blpl3 , 
Blp22,Blp23,Blp33, B2pll ,B2pl2,B2pl3,B2p22,B2p23,B2p33, 
y 1 , y2 , y3 , xl , x2 , x3 , z , sO , si , 
Degrees => {6:2, 12:1, 3:2, 3:1,2,3,4}] 

— adjointMatrix : input: M 3x3 matrix, 

output: classical adjoint of M (which is 3x3) 
adjointMatrix = (M) -> transpose matrix { 

{det submatrix(M,{l,2>,{l,2>) , -det submatrix(M, {1 , 2} , {0 , 2>) , 

det submatrix(M,{l,2},{0,l})}, 
{-det submatrix(M,{0,2},{l,2>) , det submatrix(M, {0 , 2} , {0 , 2}) , 
-det submatrix(M,{0,2>,{0,l»}, 



26 



{det submatrix(M,{0,l},{l,2}) ,-det submatrix (M, {0, 1}, {0,2}) , 
det submatrix(M,{0,l>,{0,l>)» 

— doubleAdjointMatrix: input: Ml, M2 are 3x3 matrices 

output : 3x3 matrix 
doubleAdjointMatrix = (M1,M2) -> ( 



Mill 


= Ml 


.0. 


-0; 


M112 


= Ml. 


.1 


-0; 


M113 


= Ml. 


.2. 


.0 


M121 


= Ml 


.0. 


.1; 


M122 


= Ml. 


.1 


.1; 


M123 


= Ml. 


.2. 


.1 


M131 


= Ml 


.0. 


-2; 


M132 


= Ml. 


.1 


-2; 


M133 


= Ml. 


.2. 


.2 


M211 


= M2 


.0. 


-0; 


M212 


= M2. 


.1 


-0; 


M213 


= M2. 


.2. 


.0 


M221 


= M2 


.0. 


.1; 


M222 


= M2. 


.1 


.1; 


M223 


= M2. 


.2. 


.1 


M231 


= M2 


.0. 


.2; 


M232 


= M2. 


.1 


_2; 


M233 


= M2. 


.2. 


.2 



transpose matrix { 

{M122*M233-M123*M232+ M222*M133-M223*M132 , 
-(M121*M233-M123*M231+ M221*M133-M223*M131 ), 
M121*M232-M122*M231+ M221*M132-M222*M131> , 

{-(M112*M233-M113*M232+M212*M133-M213*M132) , 
M111*M233-M113*M231+ M211*M133-M213*M131 , 
-(M111*M232-M112*M231+ M211*M132-M212*M131) } , 

{ M112*M223-M113*M222+ M212*M123-M213*M122 , 
-(M111*M223-M113*M221+ M211*M123-M213*M121) , 
M111*M222-M112*M221+M211*M122-M212*M121») 



— wedge2: input: single 3x2 matrices u, 

output: single 3x1 matrix \wedge~2 (u) 

wedge2 = (u) -> 

matrix {{det (submatrix (u, {1 , 2} , {0 , 1>) ) > , 
{-det (submatrix (u, {0 , 2} , {0 , 1}) ) } , 
{det (submatrix (u , {0 , 1} , {0 , 1>) ) » 

— mixedwedge: input two 3x2 matrices ul and u2, output 

single 3x1 matrix ul \wedge u2 + u2 \wedge ul 

mixedwedge = (ul,u2) -> 
matrix{ 

{det (submatrix (ul,{l}, {0,1}) || submatrix (u2, {2} , {0, 1}) )+ 
det (submatrix (u2, {1>, {0, 1>) || submatrix (ul, {2}, {0, 1>))>, 
{- (det (submatrix (ul,{0}, {0,1}) || 
submatrix (u2 , {2} ,{0,1}))+ 
det (submatrix (u2,{0}, {0,1}) II 
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submatrix (ul , {2> ,{0,1})))}, 
{det (submatrix(ul,{0},{0,l}) || 
submatrix (u2 , {1} ,{0,1}))+ 
det (submatrix(u2,{0},{0, 1}) || 
submatrix (ul , {1} , {0 , 1}) ) }} 
x = matrix {{xl,x2,x3}} 
y = matrix {{yl,y2,y3}} 

Al = matrix {{0,Alpl2,Alpl3},{-Alpl2,0,Alp23},{-Alpl3,-Alp23,0}} 
A2 = matrix {{0,A2pl2,A2pl3},{-A2pl2,0,A2p23},{-A2pl3,-A2p23,0}} 
Bl = matrix {{Blpll,Blpl2,Blpl3},{Blpl2,Blp22,Blp23}, 

{Blpl3,Blp23,Blp33}} 
B2 = matrix {{B2pll,B2pl2,B2pl3},{B2pl2,B2p22,B2p23}, 

{B2pl3,B2p23,B2p33}} 

— A1,A2 are generic 3x3 skew, B1,B2 are generic 3x3 symmetric 

ID = minors (2, matrix {{ yl,y2,y3, z*xl,z*x2, z*x3}, 

{ xl,x2,x3, yl,y2, y3}}) 
adjBl = adjointMatrix (Bl) ; 

adjBlpll = adjBl_0_0; adjBlp22 = adjBl_l_l; adjBlp33 = adjBl_2_2; 
adjBlpl2 = adjBl_l_0; adjBlpl3 = adjBl_2_0; adjBlp23 = adjBl_2_l; 
adjB2 = adjointMatrix (B2) ; 

adjB2pll = adjB2_0_0; adjB2p22 = adjB2_l_l; adjB2p33 = adjB2_2_2; 
adjB2pl2 = adjB2_l_0; adjB2pl3 = adjB2_2_0; adjB2p23 = adjB2_2_l; 

a = det matrix {{xl,x2,x3},{Alp23,-Alpl3,Alpl2}, 

{A2p23 , -A2pl3 , A2pl2}} ; 
invof2 = substitute ((1_S/2_S), S) ; 

b = (invof2 * (y* Bl* adjB2*Bl*transpose(y)+ 
y* B2* adjBl*B2*transpose(y)+ 
2*x*Al* adjB2* (transpose Al)*transpose(x)+ 
2*x*A2* adjBl* (transpose A2) *transpose (x) - 
2*x*Al*doubleAdjointMatrix (B1,B2)* 

transpose (A2) *transpose (x) - 
z*x*Bl*adjB2*Bl*transpose(x) - 
z*x*B2*adjBl*B2*transpose (x) + 
4*y*Bl*adj B2*transpose (Al) *transpose (x) + 
4*y*B2*adj Bl*transpose (A2) *transpose (x) - 
(Alpl2*(xl*y2-x2*yl)+Alpl3*(xl*y3-x3*yl)+ 
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Alp23* (x2*y3-x3*y2) ) * 
(Blpll* adjB2pll+Blp22* adjB2p22+ Blp33*adjB2p33 + 
2*(Blpl2* adjB2pl2+Blpl3* adjB2pl3+ Blp23*adjB2p23) ) - 
(A2pl2*(xl*y2-x2*yl)+A2pl3*(xl*y3-x3*yl)+ 

A2p23* (x2*y3-x3*y2) ) * 
(B2pll* adjBlpll+B2p22* adjBlp22+ B2p33*adjBlp33 + 
2*(B2pl2* adjBlpl2+B2pl3* adjBlpl3+B2p23*adjBlp23) ) ) ) _0_0 ; 

ul = (-A1* (transpose x) + Bl * (transpose y) ) | 

(-A2* (transpose x) + B2 * (transpose y) ) 
u2 = (-B1* (transpose x) ) | (-B2* (transpose x) ) 
sigmal = wedge2 (ul) + z * wedge2 (u2) 
sigma2 = mixedwedge (ul,u2) 

f 1 = (x*Al*transpose (y) +y*Bl*transpose (y) -z*x*Bl*transpose (x) ) _0_0 ; 
f 2 = (x*A2*transpose (y) +y*B2*transpose (y) -z*x*B2*transpose (x) ) _0_0 ; 

11 = z*xl*sO+yl*sl+sigmal_0_0; 

12 = z*x2*s0+y2*sl+sigmal_0_l; 

13 = z*x3*s0+y3*sl+sigmal_0_2; 

14 = yl*sO + xl*sl+ sigma2_0_0; 

15 = y2*s0 + x2*sl+ sigma2_0_l; 

16 = y3*s0 + x3*sl+ sigma2_0_2; 
q = sl~2-z*s0~2-a*s0+b; 

IY =ideal (f l,f2,ll,12,13,14, 15,16,q); 

The following is an example of substitution to a nongeneric case 
subslist = { xl => xl,x2 => 0,x3 => x2, 
yl => yl,y2 => y2,y3 => 0, z => 5, 
Alpl2 => l,Alpl3=>2 ,Alp23=>0, 
A2pl2 => l,A2pl3=> 3 ,A2p23=> 1, 
Blpll =>1, Blpl2 => 5, Blpl3=> 0, 

Blp22 =>1, Blp23=>l ,Blp33=> 1, 
B2pll=> -l,B2pl2 => 2,B2pl3=> 1, 
B2p22=>l,B2p23=>7 ,B2p33=>l>; 
specif icIYbef ore = substitute ( IY, subslist) ; 
specifics = kk [ yl,y2,y3,xl,x2,x3,s0,sl, Degrees => {8:1}] 
specificIY = sub (specif icIYbef ore , specif icS) 
isHomogeneous specificIY — true 
codim specificIY 
betti res specificIY 
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